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Models of medium-induced radiative parton energy loss account for the strong suppression of 
high-pT hadron spectra in ^snn = 200 GeV Au-Au collisions at RHIC in terms of a single "jet 
quenching parameter" q. The available suite of jet quenching measurements make q one of the 
experimentally best constrained properties of the hot fluid produced in RHIC collisions. We observe 
that q can be given a model-independent, nonperturbative, quantum field theoretic definition in 
terms of the short-distance behavior of a particular light-like Wilson loop. We then use the AdS/CFT 
correspondence to obtain a strong-coupling calculation of q in hot M = 4 supersymmetric QCD, 
finding ijsym = 26.69 V qsym-ZVc T 3 in the limit in which both N c and AnasYMN c are large. We thus 
learn that at strong coupling q is not proportional to the entropy density s, or to some "number 
density of scatterers" since, unlike the number of degrees of freedom, q does not grow like N% . 



Ultrarelativistic nucleus-nucleus collisions are studied 
at RHIC and at the LHC to determine the properties 
of QCD matter at extreme energy density and temper- 
ature Q, ■ If we could do the gedanken experiment of 
deep inelastic scattering on the hot fluid produced in a 
heavy ion collision, we could learn a lot. Even though 
the short lifetime of the transient dense state precludes 
the use of such external probes, a conceptually similar 
method is available at RHIC and LHC energies. This 
method is based upon internally generated probes: ener- 
getic partons produced in rare high transverse momen- 
tum elementary interactions in the initial stage of the col- 
lision, which then interact strongly with the hot, dense 
fluid produced in the collision as they plough through 
it 0- The characterization of the resulting medium- 
induced modification of high-p^ parton fragmentation 
("jet quenching") and its connection to properties of the 
hot, dense matter that is the object of study have become 
one of the most active areas of research stimulated by 
RHIC data Q. Models which supplement the standard 
perturbativc QCD formalism for high-py hadron produc- 
tion with medium-induced parton energy loss successfully 
account for the strong (up to a factor ~ 5) suppression of 
hadronic spectra in -y/ijvw = 200 GeV Au-Au collisions 
at RHIC, its dependence on centrality and orientation 
with respect to the reaction plane, and the correspond- 
ing reduction of back-to-back hadron correlations 0, 0] . 
These models typically involve one medium-sensitive "jet 
quenching parameter" denoted q. This parameter is usu- 
ally defined only perturbatively, and is often thought of as 
proportional to 1/(A|>Amfp)j with Ad the Debye screen- 
ing length and Amfp some perturbatively defined trans- 
port mean- free path 0. In the present paper, we address 
the question of how q can be defined and calculated from 
first principles in nonperturbative quantum field theory, 
without assuming the existence of quasiparticles with a 
well-defined mean-free path. 

There are many indications from data at RHIC and 



from calculations of lattice-discretized QCD thermody- 
namics that the quark-gluon plasma at temperatures 
not far above the crossover from the hadronic phase is 
strongly coupled. (For example, lattice calculations show 
that J/^S mesons remain bound Q; for example, there is 
qualitative agreement between the degree of azimuthal 
anisotropy in collisions with nonzero impact parameter 
seen in RHIC data and in calculations assuming ideal, 
zero- viscosity, hydrodynamics 

ESSES-) Lattice QCD 
is the prime example of a rigorous calculational method 
applicable in a hot, strongly coupled, gauge theory. Be- 
cause it is formulated in Euclidean space, it is well-suited 
to calculating static thermodynamic quantities and less 
well-suited to calculating transport coefficients, or dy- 
namical processes of any sort. It therefore cannot ad- 
dress parton energy loss itself, and cannot be used to 
calculate q, the property of the medium that parton en- 
ergy loss "measures" . Complementary nonperturbative 
techniques are thus desirable. One such technique is 
the AdS/CFT correspondence, which maps nonpertur- 
bative problems in certain hot strongly coupled gauge 
theories onto calculable problems in a dual gravity the- 
ory This method has been used to calculate the 
shear viscosity in several supersymmetric gauge theo- 
ries [H[l3 J9, HE El > as well as for certain diffusion 
constants lflj and thermal spectral functions E3. (See 
Refs. [13 for work towards a dual description of dynam- 
ics in heavy ion collisions themselves.) The best-studied 
example is the calculation of the shear viscosity 77, where 
the dimensionlcss ratio 77/ s, with s the entropy density, 
takes on the value l/Aix in the large number of colors 
(N c ), large 't Hooft coupling (A = 9ym.N c ) limit of any 
gauge theory that admits a holographically dual super- 
gravity description, making this result "universal" 
Furthermore, the leading 1/A corrections are known 
The results thus obtained in supersymmetric Yang-Mills 
theories have been argued to be relevant to the QCD 
matter produced at RHIC 
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In this paper, we shall demonstrate that the problem of 
calculating the jet quenching parameter also lends itself 
to a simple reformulation in theories with a gravity dual. 
After formulating a nonperturbative definition of q, we 
shall calculate q in J\f = 4 supersymmetric Yang-Mills 
theory, in the large-A c and large-A limit. We shall close 
with a comparison to RHIC data, but at present we make 
no conjecture for the "universality" of our result. 

We begin with the gedanken experiment of DIS on hot 
matter in thermal equilibrium. A virtual photon 7*, 
emitted by the electron, interacts with this matter. In 
the rest frame of the thermal matter, the incident virtual 
photon can be viewed as a superposition of its hadronic 
Fock states |y*) = \q(x)q(y)) + O(g). The lowest lying 
Fock state is a color singlet quark-antiquark dipolc whose 
transverse size in a configuration space picture is propor- 
tional to the inverse of the virtuality of the scattering 
process in momentum space. In the high-energy limit, 
the scattering S'-matrix can be written in terms of eikonal 
phase factors W, which account for the precession of the 
partonic projectile in the color field of the medium [2lj . 
For instance, for an incident quark of initial color a and 
transverse position x, 5|a(x)) = W F g(x) |/3(x)). Here 



W(-x) = V exp 



dx A + (x, x 



is the Wilson line in the representation of the projectile. 
The quark a propagates at fixed transverse position x 
along the light-cone direction x~ through the medium, 
which is described by the target color field A + . Then, 
the "photoabsorption" cross section for interaction of the 
\qq) Fock state of the virtual photon with the target is 



nis 



J d 2 xd 2 yVXx-y)V*(x-y)P t ^(x,y), ( 2 ) 



where — y) tf>* (x — y) determines the probability dis- 
tribution of color singlet dipoles of transverse size |x — y|. 
The interaction probability is [2l| 



nil 
Mot 



— Tr[W F (x)W F Hy)] 



(3) 



For a color dipole with transverse size L = |x — y| which 
propagates along the light-cone through a medium of 
length L~ , the only information about the medium that 
enters the calculation of the photoabsorption cross sec- 
tion (J2J) is then encoded in the thermal expectation value 
of a closed light-like Wilson loop {W F (C)), whose con- 
tour C is a rectangle with large extension L~ in the x~ - 
direction and small extension L in a transverse direction. 
Calculating (W F (C)) gives the entire nonperturbative in- 
put to the virtual photoabsorption cross-section J5J. 

Inspection of the calculation of the energy distribution 
of the medium-induced gluon radiation, the fundamental 
quantity of interest in jet quenching calculations, reveals 



that the medium-dependence enters this calculation via 
the expectation value of a Wilson loop in the adjoint rep- 
resentation with the same contour C [2^ . The physical 
origin of (W A (C)) in this calculation is more difficult to 
picture than was the case in the DIS example. Here, the 
Wilson loop arises in a configuration space formulation 
of the gluon radiation cross-section from the combina- 
tion of the gluon emitted at transverse position x in the 
radiation amplitude, with the gluon at position y in the 
complex conjugate amplitude. The difference L = |x — y| 
is conjugate to the transverse momentum fcj^ of the emit- 
ted gluon. (See Ref. [i^ for details.) In jet quenching 
calculations, one frequently uses the so-called dipole ap- 
proximation |23l |. valid for small transverse distances L: 



(W A (C)) 



exp 



1 
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qL~ L 2 



(4) 



Here, q is the "jet quenching parameter" introduced pre- 
viously in Ref. |6j . (The factor of y/2 in the denominator 
arises because the light-cone distance L~ is larger than 
the spatial distance travelled, conventionally used in , 
by that factor.) We simply observe that instead of seeing 
Q as an approximation, we can use it as a nonperturba- 
tive definition: —^=q is the coefficient of the L~ L 2 term 

4\/2 * 

in \og{W A (C)) at small L. 

Consider a high energy parton for which the dominant 
energy loss mechanism is the radiation of gluons with 
large enough k± ~ 1/L that QCD is weakly coupled at 
the scale k±. This energy loss process cannot be mod- 
elled in a conformal theory that is strongly coupled at all 
scales. However, the strongly coupled physics at scales 
proportional to T that characterizes the medium rather 
than the probe does enter the energy loss calculation, 
only via the quantity q. This, and the fact that time- 
averaged values of q have been determined in comparison 
with data from RHIC in Refs. 0,O|, motivates the cal- 
culation of q in strongly coupled gauge theory plasmas. 

We have calculated the thermal expectation value of 
the light-like Wilson loop 10} for Af = 4 super Yang-Mills 
theory. This is a conformally invariant theory with two 
parameters: the rank of gauge group N c and the 't Hooft 
coupling A = gy M N c . Its on-shell field content includes 
eight bosonic and eight fermionic degrees of freedom, all 
in the color adjoint representation. We begin by evalu- 
ating the li ght- like Wilson loop in the fundamental rep- 
resentation |2y|. According to the AdS/CFT correspon- 
dence j m the large- A^ c and large-A limits the thermal 
expectation value (W F (C)} can be calculated using the 
metric for a 5-dimensional curved space-time describing 
a black hole in anti-deSitter (AdS) space j^. Calling 
the fifth dimension r, the black hole horizon is at some 
r = 7'o and the (3 + l)-dimensional conformally invariant 
field theory itself "lives" at r — > 00. The prescription for 
evaluating (W F (C)) is that we must find the extremal ac- 
tion surface in the five-dimensional AdS spacetimc whose 
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boundary at r — > oo is the contour C in Minkowski space 
R 3 ' 1 . (W F (C)) is then given by 



(W F (C))=exphS(C)] , 



(5) 



with S the action of the extremal surface |24| , subject to 
a suitable subtraction that we discuss below. 

Using light-cone coordinates x M = (r, a^, x 2 , a; 3 ), the 
AdS black hole metric is given by ^lj 



cfe 2 



^ + ofe + d.T + ^ (dxj + dx 2 ) 

+\ (i~f) ((dx + ) 2 + (dx-) 2 )+jdr 2 
G^dx^dx" , (6) 



where / = 1 1 — -fij with i? the curvature radius of 
the AdS space. Here, the temperature T of the Yang- 
Mills theory is given by the Hawking temperature of the 
black hole, T = and R and the string tension l/27ra' 

are related to the t'Hooft coupling by = \/~\. 

We parameterize the surface whose action S(C) is to 
be extremized by x^ = x^ (r, a) , where a a = (r, a) de- 
note the coordinates parameterizing the worldsheet. The 
Nambu-Goto action for the string worldsheet is given by 



S = — !— / dadr y/ det g afl 
Zira J 



(7) 



with g a p = G^vdaX^dpx" the induced metric on the 
worldsheet. This action is invariant under coordinate 
changes of er Q , and we can set r = x~ and a = X2, tak- 
ing the "short sides" of C with length L along the X2 
direction. Since L~ L, we can assume that the sur- 
face is translationally invariant along the r direction, i.e. 
x n _ x ft ^ T ^ _ x n ^ _ p nc Wilson i 00 p ij es cons tant 

and at constant x + , so £3(0") = const, x + (a) = const. 
For the bulk coordinate r, we implement the require- 
ment that the world sheet has C as its boundary by 
imposing r (±|f) = 00, which preserves the symmetry 
r(a) = r(—a). The action now takes the form 



S 




2-kol'R 2 



(8) 



where r' = d a r. The equation of motion for r(a) is then 



1 R 2 



(9) 



with 7 an integration constant. Eq. © has two solutions. 
One has 7 = and hence r' = 0, meaning r(a) = 00 for 
all a: the surface stays at infinity. This solution is not 
of interest. The other solution has 7 > 0. It "descends" 
from r(±4) = 00 and has a turning point where r' = 
which, by symmetry, must occur at a = 0. From @, the 
turning point must occur where / = 0, implying r' = 



at the horizon r = r$. The surface descends from infinity, 
skims the horizon, and returns to infinity. Note that the 
surface descends all the way to the horizon regardless 
of how small L is [27j. This is reasonable on physical 
grounds, as we expect q to describe the thermal medium 
rather than ultraviolet physics. Knowing that r = ro at 
a = 0, © can be integrated, yielding 



L 
2 



7 Jr 



dr 



aR 2 
jr 



E- do) 



where a = ~ L311 - We then find 



tts/XL-LT 2 / 4a 2 
2V5 V 1 + 7r 2 T 2 L 2 



(11) 



where we have used ro = t:R 2 T and ^7- = V% 

Recalling our discussion of the DIS example, we now 
note that whereas we wish to evaluate the interaction 
of a "bare" \qq) Fock state of the virtual photon with 
the medium, what we have calculated above includes the 
"self-energy" of the high energy quark and antiquark, 
moving through the medium as if in the absence of each 
other. We perform the required subtraction upon not- 
ing that in addition to the extremal surface constructed 
above, there is another trivial one given by two discon- 
nected worldshects, each of which descend from r = 00 



to r = tq at constant x%, one at X2 



the other at 



■1-2 



The total action for these two surfaces is 



.y/XL'T 



we obtain 



Subtracting this self-energy 



Si = S — So , 



(12) 



(13) 



the L-dependent interaction between the quark pair mov- 
ing through the medium. Because the quarks are moving 
at the speed of light, both So and S are finite: no ultra- 
violet subtraction like that required for static quarks is 
needed. exp[— Si] is the thermal expectation value of the 
Wilson loop J5J in the fundamental representation, and 
would be our final result were we actually interested in 
DIS. Note that Si vanishes for L — ► as it must, since the 
photoabsorption probability |J5J vanishes in that limit. 

We now evaluate the adjoint Wilson loop Q that de- 
termines q. For SU(N C ), the Wilson line in the ad- 
joint representation can be obtained using the identity 
TrW = trWtrW' — 1, where Tr and tr denote traces 
in the adjoint and fundamental representations, respec- 
tively. This implies that in the large- N c limit, the adjoint 
Wilson loop is given by cxp [— 2Si]. The jet quench- 
ing parameter q in is determined by the behavior of 
2Sj for small transverse distances, LT <C 1. We find 



2Si 



4\/2a 



V\T 3 L-L 2 + 0(T 5 L-L 4 ) 



(14) 
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and can now use © to read off our central result: 

Qsym = —VAT 3 = r , 5 \ 4 ' VXT 3 « 26.69 V«sym iv c r 3 

(15) 

in A/" = 4 supersymmetric Yang-Mills theory in the large- 
N c and large- A limits. 

Perhaps the most striking qualitative feature of our 
result is that at strong coupling q is proportional to VA, 
not to the number of degrees of freedom ~ Y 2 . This 
means that at strong coupling q cannot be thought of 
as "measuring" either the entropy density s or what is 
sometimes described as a "gluon number density" or e 3 / 4 
as had been expected Q, since both s and the energy 
density e are proportional to V 2 A°. Instead, it appears 
that q is better thought of as a measure of T 3 |28j. 

We would be remiss not to attempt a comparison to im- 
plications of RHIC data. Taking N c = 3 and «sym = \ , 
reasonable for temperatures not far above the QCD phase 
transition, we shall use A = 6ir to make estimates. From 
(PI), we find q =4.5, 10.6, 20.7 GcV 2 /fm for T = 300, 
400, 500 MeV. In a heavy ion collision, q decreases with 
time r as the hot fluid expands and cools. The time- 
averaged q which has been determined in comparison 

with RHIC data is q = jj^yi J T ° +L ^ 2 T q( t ) dr, found 
to be around 5-15 GcV 2 /fm ^, 0|- If we assume a one- 

1 /3 

dimensional Bjorken expansion with T(t) = To (^J , 
take ro = 0.5 fm, and take L~ = 2 fm, the esti- 
mated mean distance travelled in the medium by those 
hard partons which "escape" and are detected we 
find that to obtain q = 5 GeV 2 /fm from (|15J) we need 
T such that T(l fm) w 310 MeV, only slightly higher 
than expected [8j ■ Equivalently, the q we find from i|15|) 
is slightly smaller than that suggested by RHIC data. 
First, this could indicate that in going from hot M = 4 to 
the QCD quark-gluon plasma, as the number of color ad- 
joint degrees of freedom is decreased by a factor of 2/15 
(and the color fundamental quarks are added) the jet 
quenching parameter increases slightly. Second, it could 
indicate that the q extracted from data is somewhat high, 
either because the medium through which the energetic 
parton moves has a flow velocity transverse to the par- 
ton's motion or because additional sources of energy loss 
are also significant, for example collisions without gluon 
radiation or processes occurring before the medium is in 
local thermal equilibrium. Testing the second conclusion 
likely requires further experimental data; the first can be 
tested by evaluating q in other gauge theories with grav- 
ity duals. Unless our result turns out to be "universal", 
which would be interesting in its own right, such calcula- 
tions would indicate the direction in which q changes as 
the theory is made more QCD-like. Calculating the 1/A 
corrections to our result would also be informative. 
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